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ABSTRACT 

We present an analysis of the energetics and particle content of the lobes of 24 radio galaxies at 
the cores of cooling clusters. The radio lobes in these systems have created visible cavities in the 
surrounding hot, X-ray-emitting gas, which allow direct measurement of the mechanical jet power of 
radio sources over six decades of radio luminosity, independently of the radio properties themselves. 
Using these measurements, we examine the ratio between radio power and total jet power (the radiative 
efficiency). We find that jet (cavity) power increases with radio synchrotron power approximately as 
Pjet ~ -kfadio' where 0.35 < [3 < 0.70 depending on the bandpass of measurement and state of the 
source. However, the scatter about these relations caused by variations in radiative efficiency spans 
more than four orders of magnitude. A number of factors contribute to this scatter including aging, 
entrainment, variations in magnetic field strengths, and the partitioning of energy between electrons 
and non-radiating heavy particles. After accounting for variations in synchrotron break frequency 
(age), the scatter is reduced by « 50% , yielding the most accurate scaling relation available between 
the lobe bolometric radio power and the jet (cavity) power. 

We place limits on the magnetic field strengths and particle content of the radio lobes using a 
variety of X-ray constraints. We find that the lobe magnetic field strengths vary between a few to 
several tens of microgauss depending on the age and dynamical state of the lobes. If the cavities 
are maintained in pressure balance with their surroundings and are supported by internal fields and 
particles in equipartition, the ratio of energy in electrons to heavy particles (fc) must vary widely 
from approximately unity to 4000, consistent with heavy (hadronic) jets. The synchrotron ages are 
substantially smaller than the independently measured X-ray buoyancy ages. Low frequency radio 
observations provide the best tracer of AGN power integrated over the past several hundred million 
years, and may provide a useful proxy for X-ray cavities in systems where cavities cannot be directly 
detected but are likely to exist (e.g., in cooling-flow clusters). 

Subject headings: galaxies: active - galaxies: clusters: general - X-rays: galaxies - X-rays: galaxies: 
clusters - radio continuum: galaxies 



1. INTRODUCTION 



A source of heat has been invoked to explain the lack of 
gas cooling below 2 keV at the centers of galaxy clusters 
( [Peterson et al |2001 " 



pected in the classica 



Kaast 



coolin 



ra et al.|2004[) , which was ex- 
ng flow model ( Fabianp.994 1 



sively bright cluster central galaxies ( 


Bcnson et al.||2003 


Bower et al.||2006| |Croton et al.||2006 


ISijacki & Springel 


2006 


1. The most plausible candidate is the energy in- 



jected by relativistic jets originating near super massive 
black holes (SMBH) at the centers of the clusters ( Tabor 
fc Bmney||i"993l |Binney fc Taborl[T995l |Tucker fc DavTcT 



1997k. When the jets interact with intracluster medium 
M) 



(IClVl) two lobes of plasma are inflated, which contain 
relativistic electrons and magnetic fields seen in radio 
observations. Such interactions have been discussed ex- 



tensively in the literature ( 


Scheuer 


1974[ Blandford & 


|Rees||1974| |Begelman & Cio 


m 1989; 


kaiser &; Alexander 



2 



1997 Heinz et all 1998 1 



The idea that the central active galactic nuclei (AGN) 
are the source of energy has been given impetus by Chan- 
dra X-ray Observatory images showing the keV gas be- 
ing displaced by t he radio sources harbored at the cen 



positively charge d particles beside relativi stic electrons, 
such as protons dCelotti & Fabian||1993| or positrons 

998 ^ 



(Cclotti & Fabian 
([Reynolds et al.||1996| |Wardle et aT. 
|hara||2004[ |Dunn et al.||2006ap ~' 



terof the cluster ( |McNamara et aL 2000 Fabian et al. 
2000]). Some recent analysis of samples of systems with 



X-ray cavities show that the cavities are long lasting and 
that the y provide enough energy to suppress cooling in 
clusters (jBi rzan et al l2004| [Dunn fc Fabian||2004| [Dunn 1 
et al.|2 005 ; Raffcrty ct al. 2006 ) . By consi dering only the 



cavity enthalpy, Bi'rzan et al. (20041 and Rafferty et al 
( 2006 ) show that the current outburst in more than half 



of the systems in their sample of cavity systems are pow- 
erful enough to offset cooling, although the integr ated 
energy is sufficient to o ffset cooling in all systems ( |Mc-| 
Namara & Nulsen 2007 1. However, the enthalpy is only a 
lower limit for the total energy injected with a cavity. Ev- 
idence is mounting for the existence of weak but powerful 
shocks in the ho t thermal gas. Notable examples include 
MS 0735+7421 OMcNamaraet al.|2005D, Perse us ( |Fabian| 



et al.||2006t, M8TT Forma.n ct al. 2005} |2007[ ), Hydra A 
^fn 



( |lNulsen et al.|2005b[fWise et al.|2007p among many oth 
ers. Where they have been measured, shock energies are 
comparable to cavity enthalpies. In addition, the energy 
is still probably underestimated due to further adiabatic 
losses, undetected cav ities, etc., the energy of which we 



Kino fc Taka- 
The positively charged 
particle content of the jets is observationally inaccessi- 
ble and the existence of proton-electron jets or positron- 
electron jets is inferred using indirect arguments. How- 
ever, the discovery of X-ray cavities has made it possible 
to quantify the particle content, by calculating the ratio 
between the total particle ener gy and the energy in rela- 
tivistic electron s and positro ns ( Fabian et al.|2002||Dunn| 
fc Fabian|[2004[ |Dunn et al.||2005| )~ 

In this paper, we use VLA radio data at 4 different fre- 
quencies for a sample of 24 systems with X-ray cavities to 
investigate radiative efficiencies. We derive synchrotron 
break frequencies and ages. Using both X-ray and ra- 
dio data for the cavities, we place limits on their mag- 
netic field strengths and on their particle content (i.e., 
the energies in heavy-particles and in synchrotron radi- 
ating electrons and positrons) . We investigate the effects 
that differences in the radio source ages or magnetic field 
strengths and particle content have on the efficiencies. 

Using both the X-ray and radio data, we investigate 
the radio source ages by comparing synchrotron ages 
with buoyancy ages. In general, the sy nchrotron ages are 
lower limits on the radio source age flEilek||1996j[Eilek| 



2007 see McNamara & JNuisen 2007 for a review 



et al. 1997 Eilek|[i999) due to the inhomogeneiti es m 



canno t quantify (e.g., Nusser et al. 2006 Binney et al. the magnetic held, which are difficult to est imate (Siah 



Cavities provide a lower limit to the jet energy and 
are the best- available gauges of the average jet power in 
these systems; therefore, they can be used, together with 
the radio data, to estimate the radiative efficiencies and 
radio-lobe properties. The total energy in the radio lobes 
is the sum of magnetic and particle energies: 



fc Wiita|1990| |Wiita fc Gopal-Krishna|1990[) . The X^ray 



E p = E B + (l + k)E e , 



(1) 



where k is the ratio of the energy in heavy particles to 
that in the electrons (E e ) and Eb is the energy in the 
magnetic field. The energies in the magnetic field and 
electrons are calculated as: 



derived ages (e.g., buoyancy ages) assume that the bub- 
bles detach early from the driving jet and rise buoyantly 
in the ICM. By finding discrepancies between the two 
ages, we can test the assumptions that go into their de- 
termination. Lastly, on a related topic, we are now able 
to determine more precisely whether a cavity is radio 
filled or ghost. 

We adopt H = 70 km s" 1 Mpc~\ fi M = 0.3, and 
il\ = 0.7 in all calculations throughout this paper. 

2. THE SAMPLE 



Er 



EP 
8tt 



(j>V 



E e = B 3/2 i ra dCi 2 (a, V\, v 2 ), 



(2) 
(3) 

where <fi is the filling factor (the fraction of the vol- 
ume filled by relativistic particles and mag netic field), 
Ci2(a, v\, v-i) is a constant as defined by Pacholczyk 
( 1970 1 and L ra d is the total radio luminosity: For the 
first time, the cavities allow one to estimate -Etot of ra- 
dio lobes to within a factor of a few (given the uncer- 
tainties in projection and the omission of shock energy) . 
Additionally, cavities can be used to estimate the pres- 
sures (assuming they are in pressure balance with the 
surrounding ICM) and ages of the lobes (assuming that 
the cavities rise buoyantly or at the local sound speed). 
Together, these quantities break the degeneracies in the 
energy equation and allow meaningful constraints to be 
placed on the magnetic field strengths and particle con- 
tent of the lobes. 

From radio observations we know that relativistic elec- 
trons and magnetic field are present in the jets. Un- 
der the assumption of charge neutrality, there must be 



We use the Birzan et al. (20041 sample of 18 systems 
with well-defined surface brightness depressions associ- 
ated with their radio source, plus 6 new cavity sys- 
tems recently discovered by Chandra (Sersic 159/03, 
Zw 2701, Zw 3146, MACS J1423.8+2404, A1835 and 
MS0735+74). The sample consist of 22 galaxy clusters, 
one galaxy group (HCG 62), and one elliptical galaxy 
(M84). There is a broad range in redshift from 0.0035 
(M84) to 0.545 (MACS J1423. 8+2404) and a broad range 
in radio luminosities from 1.5 x 10 38 erg s _1 (HCG 62) 
to 6.9 x 10 45 erg s _1 (Cygnus A). The X-ray proper- 
ties (pressure, volume, cavity power, cavity ages) that 
we used in this work are from |Rafferty et al. ( 2006 1 . 



3. RADIO OBSERVATIONS AND DATA REDUCTION 

In order to have radio data at 327 MHz, 1.4 GHz, 4.5 
GHz and 8.5 GHz we observed 13 objects that did not 
have VLA archival data at one or more of these frequen- 
cies. The observations that we used for each frequency 
are given in the Appendix. Data were analyzed with the 

1 The VLA (Very Large Array) is a facility of the National Ra- 
dio Astronomy Observatory (NRAO). The NRAO is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Astronomical Image Processing System (AIPS). The pri- 
mary flux density calibrators and the phase calibrators 
are also listed in the Appendix. The data were Fourier- 
transformed, cleaned, restored, and self-calibrated with 
AIPS. In order to subtract confusing sources we imaged 
the entire primary beam. We list in the Appendix the 
resolution, position angle and the RMS noise of the final 
radio images for at least three different frequencies (327 
MHz, 1.4 GHz, 4.5 GHz and/or 8.5 GHz). Where avail- 
able, data at different arrays were added for better u-v 
coverage. In some cases we show images at different res- 
olutions for the same frequency: A133 at 330 MHz and 
1.4 GHz, RBS 797 at 1.4 GHz and 8.5 GHz, and A1835 
at 1.4 GHz. 

For each source we measured the total flux density us- 
ing TVSTAT in AIPS. In order to measure the flux den- 
sity for the core of the radio sources, we used JMFIT in 
AIPS to fit a gaussian fixed to the beam shape. This 
procedure was performed only for those images in which 
the lobes are clearly resolved and the core is evident. 
For these systems, the flux in the lobes (including any 
jet emission) was calculated by subtracting the core flux 
density from the total flux density. In Table [T] we list the 
measured flux densities at 327 MHz, 1.4 GHz, 4.5 GHz, 
and 8.5 GHz for the total source for all the objects in our 
sample and for the lobes for the subsample of objects for 
which the core could be accurately subtracted. We note 
that for this subsample, the core's contribution to the 
total flux is small (< 20%). The errors in the flux densi- 
ties include multiplicative error s plus RMS noise in the 
radio map (additive errors). In |Carilli et al. (|1991t the 



multiplicative errors in the flux densities range between 
2% at 4.995 GHz and 6% at 327 MHz. For simplicity, 
we chose an average of 4% for the multiplicative errors 
at all frequencies. Thi s value is consis tent with the mul- 
tiplicative error that Slee et al. (2001) used for the 1.425 
GHz flux density (4.5% to 5.0%J7^ 

4. RESULTS AND DISCUSSION 

4.1. Radio Morphologies 

The radio contour maps for all the systems listed in 
Table [T] are shown in the Appendix. Where there are 
2 images for the same frequency we list the arrays that 
were used for each image. In some systems, we do not 
see extended emission at any of the 3 or 4 frequencies 
that were observed (e.g., Zw 3146). In general, the lack 
of extended radio emission is seen at higher frequencies 
(4.5 GHz and 8.5 GHz) where for some of the systems we 
imaged only the core or the inner lobes (e.g., Hydra A; 
see Table [I] where we list the core and lobe contributions 
for all the systems). Since the array configurations were 
chosen so that the largest angular scale of the resulting 
image would be greater than the size of the X-ray cavity, 
the lack of extended radio emission at a given frequency 
is likely due to the observed frequency being above the 
break frequency of the lobes (and thus any emission is 
below the sensitivity of the observations) and not due to 
a lack of short antenna spacings. On the other hand, the 
majority of the lower- frequency observations (327 MHz 
and 1.4 GHz) show extended emission in the lobes. Ad- 
ditionally, the lobes of some systems are barely resolved, 
especially in the lower-frequency observations (e.g., A478 
at 327 MHz). 

By examining the radio morphology at different fre- 



quencies, we confirm that two of the systems in our sam- 
ple have jets or lobes oriented differently at higher fre- 
quencies (4.5 GHz and 8.5 GHz) than at lower frequencies 
(327 MHz and 1 .4 GHz). Systems with this morphology 



are Abell 2597 (Polla ck et al... 
and RBS 79 7 ( |Gitti et 

d that di 



2005 Clarke et al. 



Additionally, 



2005) 



et al.| ( |2006b[ ) demonstrated that different generations of 
X-ray cavities in the same cluster are not always observed 
at the same position angle with the respect of cluster 
core. This morphology can be explained by either a de- 
flected jet ([Pollack et al j|2 005|, by jet prece ssion ( |Dunn 
et al ||2006b n, or by both (IGitti et al j[2006[ ). As an ex- 
a mple, the S' shape jets in Abell 1795 were explained 
by Rodriguez-Martinez et al. (20061, using hydrodynam- 
ical calculations, as a point-symmetric precessing jet in- 
teracting with a plane-parallel wind. However, because 
our data were obtained to image extended lobe emission, 
they are usually of insufficient resolution to resolve the 
jets in the sample as a whole. 

A persistent problem with cavity heating in simula- 
tions is the distribution of energy, as an unchanging jet 
orientation will tend to heat the ICM anisotropically, 
forming low-density chan nels that carry the jet's ene rgy 
out of the cooling region ( Vernaleo fc Reynolds|2006" cf. 
Heinz et al.||2006 l. However, if such differently oriented 
jets are common, the cavity energy might naturally be 
spread isotropically t hroughout the atmosphere o f a clus- 
ter (for a review, see McNamara & Nulscn 20071. 



4.2. Spectra and Total Radio Luminosity 

The synchrotron spectrum emitted by an ensemble 
of electrons is well characterized by an initial power- 
law electron energy distribution N(E) cc E \ with 
the corresponding synchrotron spe ctrum I{v) oc v~ a , 
where 7 = 2a + 1 ( |Pacholczyk| 1 970' I . Synchrotron losses 
steepen the spectrum at higher frequencies by depleting 
high energy electrons first. A curvature or break occurs 
in the spectrum at the critical frequency, vq, which is 
d etermined b y synch rotron age. 

Kardashevj (1962) studied the synchrotron spectrum 
of an ensemble of electrons evol ving with time in the 
presence of a magnetic field (see Pacholczyk 1970 1 and 
showed that the steepening of the electron spectrum at 
the high frequencies is dependent on the pitch angle and 
energy of the relativistic electrons. There are three main 
models for characterizing the synchrotron spectrum: the 
KP (Kardashev-Pacholczyk) model, which assumes a sin- 
gle injection of electrons and a pitch angle for the radi- 
ating electrons that is constant with time; the JP (Jaffc- 
Perola) model, which assumes a single injection of elec- 
trons and an isotropic pitch angle distribution; and the 
CI (continuous injection) model, which assume s a mixed 
population of electrons of different ages (see Myers fc 



ack et al 
The brea 



1998) 



Spanglcr 1985; |Carilli et al.|[l99l] |Feretti et al.||1998| 

IT 



i frequency is determined by fitting these 
models to the rad io spectrum (for details of the code, see 
Carilli et al. 1991 ). The free parameters of the models are 
the injection index, a; n j , which determines the power law 
index at frequencies below the critical frequency, and the 
critical frequency itself, vq, above which the spectrum 
steepens. We found that the CI model did not provide 
acceptable fits to many of our spectra, so we do not con- 
sider it further. We fit the data with both the KP and JP 



TABLE 1 
Radio fluxes and spectral indices. 









Component 


(^327 
(Iv) 


\ J J J 


■ i ~i mi 


5 gron 




a 




A133 





060 


lobes 


3.5 ± 0.2 


0.131 ± 0.006 


















total 


3.6 ± 0.1 


0.132 ± 0.005 




0.0029 ± 0.0001 


2.26 


± 


0.04 


A262 





016 


lobes 


0.244 ± 0.015 


0.047 ± 0.004 


















total 


0.299 ± 0.012 


0.073 ± 0.003 




0.0042 ± 0.0002 


0.98 


± 


0.04 


Perseus 





018 


total 


24.5 ± 1.0 


23.2 ± 0.9 




23.9 ± 1.0 


0.03 


± 


0.03 


2A 0335+096 





035 


total 


0.21 ± 0.01 


0.037 ± 0.002 b 


0.010 ± 0.001 


0.004 ± 0.0002 


1.2 


± 


0.1 


A478 





081 




0.11 ± 0.01 


0.027 ± 0.001 




0.0054 ± 0.0002 


0.96 


± 


0.07 


MS 0735.6+7421 





216 


lobes 


0.72 ± 0.04 


0.0117 ± 0.0012 


















total 


0.80 ± 0.03 


0.021 ± 0.001 




0.0013 ± 0.0001 


2.48 


± 


0.04 


Hydra A 





055 


outer lobes 


36 ± 10 


















inner lobes 


116 ± 5 


39.1 ± 1.6 


14.8 ± 0.6 


7.2 ± 0.3 














total 


152 ± 6 


39.2 ± 1.6 


15.0 ± 0.6 


7.4 ± 0.3 


0.93 


± 


0.04 


RBS 797 





350 




0.104 ± 0.006 


0.021 ± 0.001 


0.0042 ± 0.0003 


0.0026 ± 0.0001 


1.08 


± 


0.05 


Zw 2701 





214 


total 


0.21 ± 0.01 




0.0043 ± 0.0002 


0.0029 ± 0.0001 


0.7 


± 


0.1 


Zw 3146 





291 


total 


0.028 ± 0.003 




0.00139 ± 0.00007 


0.00078 ± 0.00005 


1.04 


± 


0.14 


M84 


0.0035 


lobes 


10.9 ± 0.5 


5.38 ± 0.23 


2.07 ± 0.09 


1.26 ± 0.07 














total 


11.1 ± 0.5 


5.6 ± 0.2 


2.28 ± 0.09 


1.51 ± 0.06 


0.46 


± 


0.04 


M87 


0.0042 


lobes 


119 ± 5 


133 ± 6 


54 ±3 


44 ± 2 














total 


124 ± 5 


138 ± 6 


59 ± 2 


47 ± 2 


0.68 


± 


0.03 


Ccntaurus 





Oil 


lobes 


7.1 ± 0.7 


2.82 ± 0.16 


1.24 ±0.06 


0.77 ± 0.04 














total 


12.3 ± 0.5 


3.4 ± 0.1 


1.37 ± 0.06 


0.84 ± 0.03 


0.82 


± 


0.04 


HCG 62 





014 


total 


0.008 ± 0.002 


0.0050 ± 0.0004 




0.0015 ± 0.0001 


0.33 


± 


0.16 


A1795 





063 


total 


3.36 ± 0.14 


0.88 ± 0.04 




0.099 ± 0.004 


0.89 


± 


0.04 


A1835 





253 


total 


0.095 ± 0.007 


0.031 ± 0.001 


0.0099 ± 0.0004 


0.0074 ± 0.0003 


0.77 


± 


0.06 


MACS J1423.8 





545 


total 


0.027 ± 0.002 


0.0044 ± 0.0002 






1.22 


± 


0.06 


A2052 





035 


lobes 


27.3 ± 1.4 


5.3 ± 0.3 


0.39 ± 0.04 
















total 


30.9 ± 1.2 


5.7 ± 0.2 


0.72 ± 0.03 


0.264 ± 0.012 


1.12 


± 


0.04 


MKW 3S 





045 


lobes 


4.71 ± 0.18 


0.08 ± 0.04 


















total 


4.73 ± 0.18 


0.090 ± 0.004 


0.0025 ± 0.0001 


0.0021 ± 0.0001 


2.66 


± 


0.04 


A2199 





030 


lobes 


24 ± 1 


1.64 ± 0.07 


0.16 ± 0.02 


0.07 ± 0.01 














total 


24 ± 1 


3.68 ± 0.12 b 


0.31 ± 0.01 


0.152 ± 0.006 


1.28 


± 


0.04 


Cygnus A 





056 


lobes 




1220 ± 68 


370 ± 23 


123 ± 10 












total 


4375 ± 194 


1450 ± 60 


475 ± 20 


180 ± 10 


0.72 


± 


0.04 


Sersic 159/03 





058 


total 


1.53 ± 0.06 


0.22 ± 0.01 


0.056 ± 0.002 


0.024 ± 0.001 


1.31 


± 


0.04 


A2597 





085 


total 


8.3 ± 0.3 


1.86 ± 0.07 


0.37 ± 0.02 


0.118 ± 0.005 


1.04 


± 


0.04 


A4059 





048 


lobes 


9.0 ± 0.4 


0.83 ± 0.04 


0.063 ± 0.003 


0.016 ± 0.001 














total 


9.93 ± 0.40 


1.28 ± 0.04 b 


0.071 ± 0.003 


0.020 ± 0.001 


1.38 


± 


0.04 



Note. — When we have multiple final images of different resolutions (sec Table m), the lower resolution image was used to measur e 
the flux, with the exceptions of A133 and Hydra A. Hydra A: 74 MHz and 327 MHz radio maps are fro m |Lane et al.| ( 2004); 
g-total _ 441 _j_ 20 Jy; ^mncriobos _ 2ig ± 1Q Jy . ^outcriobcs _ 2 2l ± 30 Jy. A2597: 327 MHz radio map is from |Clarke et al.|j2005( . 
a The spectral index, a was calculated using the fluxes at the two lowest available frequencies assuming a power-law spectrum oc 



The fluxes arc from literature: sec Table [2] for references. 

models at a given injection index and selected the break 
frequency that corresponded to the best fit of the model 
to the data. In general, the KP and JP models give sim- 
ilar fits of the lobes' spectra, but in some cases the KP 
model fits better the higher-frequency points. Therefore, 
for simplicity, we use the break frequencies derived with 
the KP model. 

Table [l] gives the radio fluxes for the lobes and to- 
tal source and the spectral index for the total source, 
denoted by a (calculated between the two lowest fre- 
quencies and used later in Section |4.7[ ) . Table [2] lists 
the break frequencies for the lobes and total source that 
we derived by fitting the spectra (shown in Figure nl 
with a KP model, the derived injection index, and the 
corresponding chi-squared statistic. For our sample, the 
injection index varies from 0.5 to 1.3. However, the pre- 
ferred value for the injection index in the literature is 0.5 
or sometimes 0.75 (Komissarov & Gubanov 19941. Our 
higher values of the injection index for some objects sug- 
gest that our data are sensitive to a different region of 
the asymptotic spectrum than those of Slee et al. ( 2001[ ) 
and Komissarov & Gubanov (19941. As a result, these 
sources may have a much lower break freque ncy than we 



measure, consistent with previous fi ndings (Komissarov 
k, Gubanov] |1994| |Slee et al.||2001| ). For a number of 
systems, we were not able to obtain lobe fluxes at three 
frequencies (see Table [lj. In these cases, to derive the 
break freq uency , we use the approximation of M yers fc] 



in the last column of Table |2] 

We calculated the bolometric radio luminosity for the 
total source and the lobes by integrating under the best- 
fit KP model (given by the ai n j and vq listed in Table 
2| between the rest-frame frequencies of 10 MHz and 
TOOOO MHz. Errors in the resulting luminosities include 
errors on vq and the model normalization. In cases where 
vq is an upper or lower limit, we calculated the errors 
on the luminosity assuming that vq could be as low as 
100 MHz or as high as 10000 MHz. In all cases, we 
symmetrized the errors by adding the upper and lower 
errors in quadrature. 

We note that the lower cutoff frequency of 10 MHz 
that we adopt, while commonly used in the literature, 
is arbitrary and results in an arbitrary energy cutoff for 
the electrons. Reducing the low-energy cutoff will tend 
to increase the electron pressure which in turn will re- 
duce k. This effect would be greatest for the sources 
with steepest spectra at low frequencies (largest cti n j)- 
However, while lowering the cutoff frequency for these 
steep spectrum sources might reduce the large k values, 
our data do not extend to such low frequencies, and we 
cannot therefore choose a cutoff frequency based on our 
observations. For simplicity and consistency with much 
of the literature, we adopt a cutoff of 10 MHz and note 
that further investigation of the low-energy cutoff is war- 
ranted. 



Spangler (1985) using the fluxes at the frequencies given 
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Fig. 1. — The spectra of the radio lobes, roughly organized from the flattest spectrum (top left) to the steepest spectrum (bottom right). 
Each spectrum is shifted by an arbitrary factor so that all spectra can be shown in one plot. In order to have as many data points as 
possible, in some cases we complemented our data with literature values (see Table p|. The data were fit with the KP model in order to 
derive the break frequency. We show the fit with the best \ 2 value. Both the input parameters for the KP model («i n j) and the results of 
the model (uc, X 2 ) are listed in Table [5] 



4.3. Cavity Classification 

In the literature, X-ray cavities have been divided into 
two categories depending on the presence of bright 1400 
MHz radio emission in the cavities: radio-filled and radio- 
ghost cavities. Radio-filled cavities, such as the inner 



cavities in Hydra A ( |McNamara et al. 2000 1 , have 1400 
MHz or higher freq uency emission. Gh ost cavities, such 
as those in Perseus ([Fabian et al.||2000| ) and A2597 ( |Mc- 
Namara et al. 2001 1, lack significant emission at 1400 
MHz. The ghost cavities were interpreted as detached 
relics of earlier outbursts whose radio emission has faded 
(McNamara et al. 20011. However, this classification 



depends critically on the available radio images and is 
largely a qualitative one. In this section we give a quanti- 
tative classification system based on the break frequency 
of the synchrotron emission. 

We have radio data at 4 different frequencies for most 
of the systems in our sample (see the Appendix) . There- 
fore, we can now determine more precisely whether a 
cavity is radio filled or ghost. Overlays of the radio emis- 
sion (in green) on the smoothed Chandra X-ray images 
are shown in the Appendix. Since we are interested in 
the correlation between the X-ray cavities and the ra- 
dio emission at different frequencies, we display only the 
cases where we have extended radio emission. These fig- 
ures show that Hydra A (inner lobes), M84, M87, Cen- 
taurus, A2052, A2199, and Cygnus A have X-ray cavities 
filled with 1.4 GHz and higher frequency radio emission 



(4.5 GHz and/or 8.5 GHz) and would be classified as 
radio-filled cavities based on the above classification. On 
the other hand, A133, A478, and MS 0735.6+7421 were 
classified as radio-filled cavities based on the fact that 
1.4 GHz fills the cavities; however, our new data do not 
show any higher frequency emission filling the cavities in 
these systems. 

Since all the cavities have some radio emission when 
probed deeply, the traditional distinction between radio- 
ghost and radio-filled cavities is somewhat arbitrary. In 
order to be more quantitative, we propose that the break 
frequency be used to separate cavities between radio- 
filled and radio-ghost. The traditional classification is 
essentially equivalent to saying that a radio-filled cavity 
is one that has a break frequency higher than 1.4 GHz, 
since the break frequency defines the point above which 
the spectrum falls steeply. With this idea in mind, we 
classify systems with lobe break frequencies (yc) below 
~1.4 GHz as ghost cavity systems and those with break 
frequencies above 1 .4 GHz as radio-filled systems (see Ta- 
ble 2J. However, note that many systems are classified 
as ghost because of the lack of a lobe break frequency es- 
timate (e.g., Zw 2701, Zw 3146, etc.). Also, we note that 
1 .4 GHz is an arbitrary frequency for separating the sys- 
tems between radio- filled and ghost, since Table [2] shows 
a continuous range of break frequencies. However, as we 
show in Section |4.4| there is a clear separation between 
the synchrotron ages of the radio-filled and radio-ghost 
systems, consistent with the interpretation that the ghost 
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324, 408 (10), 4860, 8460 
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total 
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References. — (1) |Andernach et al.Hl981>; (2) |Bennett| l|1962t; (3) |Burbidge fc CrowneHl979); (4) | Condon et al.| jl998|l; (5)|Edge et al.|11959 
(6)|Gower et al.| l|1967t; (7)|Cregory & Condon| Q199H ; (8) lAales et al.| 1199H; (9) |Haynes et al.|jl975{; (10)|Large et al.|]199lt; (11) |Pauliny- l'ot 



fll967TT(7)|Uregory & Condon| | |199H; (8)|Hales et al.|1199U: (9) | Haynes etafTi|1975t; ( 10)|Large et al-H1991>; (ll)|i 
I) IPtlkington fc ijcottj jj^65j; (lit) |tjadlerHi^4t; ( TV) |}jlee| (1995) ; (15) |Vollmer eTat] l |2005[ ; (lti) |Wall jit Peacock| 



et al.HiMi); jf j 

Waldram et al.| |l996p ; (18) |Wright fc Utrupcek| fT990p ; (19) TWright et al.|]199 



TTT> 



The injection index used in the KP model. 

b The reduced chi-squarcd statistic given by the KP model. 

c The break frequency corresponding to the best-fit KP model. Those values marked with an asterisk were computed using the relation of Myers & 
|Spangler| l |1985| , assuming Qi n j — 0.5. 

The linear expansion factor for a volume of plasma. 
c The adiabatic-loss-corrcctcd break frequency. 

f The frequencies in MHz for the radio data used to derive the break frequency (i^c); the number in parentheses are the references for the literature 
values of the radio fluxes. 



cavities are relics of an earlier outburst. 

However, the situation has become complicated with 
the discovery of tunnels in the X-ray emission filled with 



the lo bes is fading away (as su ggested by Heinz et al 
2002}. A2597 is interpreted by IClarke et al 



MS 



low frequen cy radio emission (Clarke et al. 2005 



et al.|2007| . Such structure is possibly indicative of con 



Wise 



tinuous feeding of th e lobes by the cent ral source, such similar situation ( Wise et al 



as in the models of Dichl et al. (2008). Additionally. 



an intermediate case based on the finding of an X-ray 
tunnel which connects the western ghost bubble to the 
core of the source. The outer lobes in Hydra A are in a 

2007| . 



some of the ghost systems show inner jets or lobes at 
higher frequencies (Hydra A - outer lobes, RBS 797, Ser- 
sic 159/03, A2597, and A4059). These radio sources can 
be either restarting or fading away. For RBS 797, based 
on the discovery that the inner jets seen at higher fre- 
quencies have a different orientation than the outer lobes 
seen at lower frequencies, the interpretation is that this 



For Sersic 159/03, 

the 10 ks X-ray observation was not deep enough to im- 
age the cavities well. The 327 MHz radio map shows 
extended emission, but o riented different l y tha n the X- 



ray cavities identified in Rafferty et al. (2006). Much 



radio source has restarted (see also Gitti et al. 20061. 



However, for A4059 the inner jets have the same orienta- 
tion as the outer ones; therefore, this radio source may be 
in an intermediate stage where the radio emission from 



deeper X-ray observations are required in order to clar- 
ify this discrepancy. However, the radio source seems to 
be either restarted or fading away; as a result the X-ray 
cavities are classified as ghost. 

4.4. Synchrotron Ages 

Under the assumption that the relativistic electrons in- 
side the radio lobes lose their energy only by synchrotron 
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emission and inverse Compton scattering, for the KP 
model the synchrotron break (critical frequency) and the 
synchrot ron age of a distribution of electrons are rela ted 
through ( |Slee et al.|[200l| |Feretti fc Giovannini||2007[ ): 

rjl/2 

t fivn = 1060 — tttttMI + z)]- 1/2 Myr, (4) 



B 2 



2 B 2 
3 m 



where the magnetic field B is in [iG and the break fre- 
quency vq is in GHz. The magnetic field equivalent of 
the 3K background radiation, B m , accounts for inverse- 
Compton losses, with B m = 3.25[1 + z] 2 fiG. 

The largest uncertainty in the syn chrotron age calcula- 
tion is the magnetic field strength. 



Siah & Wiita ( 1990 1 



and Wiita & Gopal- Krishna (1990) discussed the diffi- 



culty of estimating the strength of the magnetic field 
and how inhomogeneities in the ma gnetic field can af- 
fect the spectral age estimates. Also, Tribble (19931, by 
studying the effects of a random magnetic field on the 
spectral age calculation, showed that the magnetic field 
rather than the pitch angle distribution is responsible for 
the diffi culty of esti mating the spectral age (JP or KP 



model) . Eilek (I1996J showed that inhomogeneities in the 
magnetic field can have important consequences for the 
radio sources, e.g., produ cing a deceptiv ely short radia- 
tive age (for a review, see Rudnick|2002 1 . 

We computed the synchrotron ages of the lobes using 
equation [4] and the minimum energy estimate of the mag- 
netic field strength in the lobes, B eq , or the magneti c field 



4.6.1 



required for pressure equilibrium, B p (see Section 
for a description of how the magnetic field strengths were 
calculated). In Table|3]we list the equipartition magnetic 
field strengths, the magnetic field strengths required for 
pressure equilibrium (for k — and 0=1), and the cor- 
responding synchrotron ages (t^ n , if yn ). We note that 
there is a clear separation in age between those cavities 
classified as radio filled (t syn < 0.1 x 10 7 yr) and those 
classified as radio ghost (t syn > 2 x 10 7 yr) , a result which 
tends to support our classification scheme. 

"Spectral ages" computed using equation El (t^) are 
estimates of the time elapsed since acceleration of the 
particles, under the assumption that synchrotron radia- 
tion is the only mechanism by which the electrons lose 
energy. However, in addition to synchrotron losses, adia- 
batic expansion can also deplete particle a nd field energy 
in a radio lobe ( Scheuer & Williams 1968 ) , reducing the 
break frequency and the magnetic field, so making the 

exceed the true age. 



measured spectral age, t^ D 



Blundell & Alexander (19941 computed a formula for 
break frequency versus age, which takes into account 
both synchrotron losses and adiabatic-expansion losses: 
vq = i/™ r /(A), where vq is the measured break fre- 
quency in the presence of adiabatic-expansion losses, and 
/(A) = 16A 6 /[(A + 1) 2 (A 2 + 1) 2 ], where A is the linear- 
expansion factor for a volume of plasma. We calculated 
A as the ratio between the final radius of the bubble at 
the end of the adiabatic expansion (r\) and the initial 
radius of the bubble at the start of the adiabatic ex- 
pansion (r ). According to the adiabatic-expansion law: 

r\ = tq(pi/pq)~^ , where 7, the ratio of specific heats, is 
assumed to be 4/3, and p\ and po are the gas pressures at 
the location of the X-ray cavity and at the center of the 
cluster respectively ( derived from X-ray observations, see 
Rafferty et alp)06] ). We assume that the bubble starts 



at the cluster center and calculate r$ from this equation. 
Table [3] lists the adiabatic-loss-corrected ages (tsy n:Cor ), 
and the factors involved in its derivation (A, /(A) and 
v c ° r ) are listed in Table [2 Our assumption that the par- 
ticles are all injected at the central pressure and expand 
to the current lobe pressure probably overestimates the 
adiabatic correction, which is generally modest. 



4.5. Comparing Age Estimates 

In Figure [2] we compare the adiabatic-loss-corrected 
synchrotron ages of the radio l obes to the X-ray-der ived 
ages of the cavities (taken from Rafferty et al.|[2006 and 
averaged over each set of cavities). There is no strong 
trend in this figure between the two age estimates, im- 
plying that they are in general only weakly coupled. Ad- 
ditionally, the clear separation between radio-filled and 
radio-ghost cavities discussed in the previous section is 
evident in this figure. 

When comparing the two estimates of age, one can 
assume either that the cavity ages are the true ages or 
that the synchrotron ages are the true ages. If we as- 
sume that the cavity ages represent the true ages, we 
find that the synchrotron ages systematically underesti- 
mate the ages of the radio sources by a mean factor of 
~ 13 with a median value ~ 1.5 when the buoyancy es- 
timate is used. When the sound-speed age is used, the 
situation improves somewhat, although the radio- filled 
cavities still have synchrotron ages that are much lower 
than their cavity ages. We can also constrain the ages in 
the radio-filled systems from the weak shocks seen in X- 
ray observations b y fitting models to the surface bright- 



ness profiles (e.g., Nulsen et al. 2005b|a l. For example 
the shock age for Cygnus A based on the surface bright- 
ness profile from Smith et al. (20021 is around 1.7 xlO 7 
yr. As a result, the buoyancy ages m radio-filled systems 
are probably several times too large, but the synchrotron 
ages are an order of magnitude too sm all. In genera l, our 
results agree with the conclusions of Eilek (11999), who 
finds that the synchrotron ages are generally lower limits 
to the radio source ages, and support the idea that parti- 
cle injection is still occurring in the radio- filled systems. 

If, however, we assume that the synchrotron age rep- 
resents the true dynamical age of the cavity system, Fig- 
ure [2] implies that many of the cavities, particularly the 
radio-filled ones, are being strongly driven, or that they 
form at radii close to their observed radii. If we assume 
that they are driven, Figure [2] implies that they are mov- 
ing ~ 4 to ~ 20 times faster than the buoyancy age sug- 
gests, so in most cases substantially supersonically. As a 
result, in these systems we would expect to see sh ocks. 
The discovery of mild sh ocks in some systems (see Mc- 
Namara & Nulsen |2007[ ) fits with the above scenario, 
and we may expect to see shocks in all of the other sys- 
tems for which the synchrotron age is much less than the 
buoyancy age of the cavity. 

We can test in detail whether the synchrotron ages 
are indicative of the dynamical ages by calculating the 
implied Mach numbers, M = R/(t^ ncoT c s ), where R is 
the distance from the cluster center to the cavity cen- 
ter and c s is the velocity of sound in the pre-shock gas, 
both taken from Rafferty et al. (20061. A number of ob- 
jects have high Mach numbers, between 3.5 and 8.0 for 
M84, M87, and Centaurus, 14 for the inner lobes of Hy- 
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TABLE 3 

Lobe properties derived using the cavity sizes from X-ray observations. 









B p c 


+eq d 

syn 




tj J 


B buoy s 


1 ~T~ &cq 


1 ~f"fe buoy 




Name 


Class a 


(/xG) 


(/iG) 


(10 yr) 


(10 7 " yr) 


(10^ yr) 


(/iG) 






(/iG) 


A133 


RG 


8.45 ± 0.14 


54 ± 1 


> 5.8 


> 4.8 


> 0.43 


13.8 ± 0.2 


340 ± 20 


180 ± 10 


45 ± 1 


A262 


EG 


4.3 ± 0.1 


44 ± 6 


10.3 ± 0.1 


7.7 ± 0.1 


0.4 ± 0.2 


24.0 ± 0.5 


1800 ± 800 


2200 ± 400 


37 ± 5 


MS 0735+7421 


RG 


4.7 ± 0.1 


42 ± 3 


9.7 ± 0.1 


5.0 ± 0.7 


0.6 ± 0.1 


4.2 ± 0.1 


1100 ± 300 


150 ± 20 


35 ± 3 


Hydra A-outcr J 


RG 


3.2 ± 0.1 


29 


25 ± 1 


10.2 ± 0.3 


1.6 


7.5 


1300 


400 ± 200 


24 ± 3 


Hydra A-inncr 


RF 


32 ± 1 


74 ± 2 


0.221 ± 0.006 


0.135 ± 0.004 


0.06 ± 0.01 


6.6 ± 0.4 


13 ± 1 


NS k 


69 ± 3 


M84 


RF 


35 ± 4 


NS k 


< 0.2 


< 0.1 




25.9 ± 0.5 


NS k 


NS k 


40 ± 20 


M87 


RF 


80 ± 10 


NS k 


< 0.06 


< 0.05 




25 ± 1 


NS k 


NS k 


90 ± 10 


Ccntaurus 


RF 


33 ± 4 


67 ± 3 


< 0.21 


< 0.15 


< 0.07 


13.8 ± 0.8 


6.8 ± 1.1 


2.6 ± 2.3 


60 ± 10 


A2052 


RG 


14.2 ± 0.1 


64 ± 1 


1.9 ± 0.1 


1.6 ± 0.9 


0.20 ± 0.01 


22 ± 1 


100 ± 10 


80 ± 20 


53 ± 2 


MKW 3S 


RG 


10.1 ± 0.2 


53 ± 3 


> 5.2 


> 4.0 


> 0.4 


7.5 ± 0.1 


180 ± 30 


39 ± 3 


45 ± 2 


A2199 


RG 


11 ± 2 


70 ± 2 


> 4.0 


> 2.6 


> 0.3 


18.2 ± 0.5 


300 ± 10 


160 ± 20 


58 ± 2 


Cygnus A 


RF 


56 ± 5 


74 ± 1 


0.15 ± 0.01 


0.077 ± 0.004 


0.10 ± 0.01 


5.0 ± 0.2 


1.8 ± 0.1 


NS k 


76 ± 6 


A4059 


RG 


13.5 ± 0.6 


46 ± 3 


3.5 ± 0.1 


2.9 ± 0.8 


0.56 ± 0.01 


17 ± 4 


41 ± 9 


33 ± 26 


40 ± 3 



a The cavity classification: radio-filled cavites (RF) have lobe break frequencies (vc, see Table |2| above 1.4 GHz; radio-ghost cavities (RG) have lobes 
with break frequencies below 1.4 GHz. 

b The equipartition magnetic field strength for k — and — 1. 

c The magnetic field strength required for pressure equilibrium for k — and = 1. 

d The synchrotron age computed using the equipartition magnetic field strength (i? eq ) and the break frequency uq. 

e The adiabatic-loss-corrcctcd synchrotron age computed using the equipartition magnetic field strength (-B cq ) and the adiabatic expansion corrected 
break frequency u%? r . 

f The synchrotron age computed using the magnetic field strength required for pressure equilibrium (B p ) and the break frequency vq. 
s The magnetic field strength obtain by assuming that the equipartition synchrotron ago (t°y n COT ) is equal to the buoyancy age. 

Assumes 0—1. NS denotes cases for which no physical solution was possible. 
1 The magnetic field strength that corresponds to the fc cq value for which the equipartit ion magnetic field equals the field required for pressure equilibrium. 
J The X-ray volume, buoyancy age, and sound speed age are from |Wise et ah] ]2007[ . 
k No physical solution 
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Fig. 2. — The synchrotron age of the lobes corrected for adiabatic-expansion losses versus the X-ray— derived buoyancy age (left) and 
sound-speed age (right) of the X-ray cavities. Triangles denote radio-filled cavities, and circles denote radio-ghost cavities (see Section |4.3| 
for details). Dashed lines denote different ratios of the synchrotron age to the cavity age. 



dra A, and 24 for Cygnus A. Such Mach numbers imply 
strong shocks that arc not seen; the Mach numbers of 



ages; see Murgia et al. 1999 1. For the ghost cavities, the 
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one may expect to see shocks whenever the synchrotron 
age is much lower than the buoyancy age. 

Therefore, we conclude that the synchrotron and buoy- 
ancy ages appear to be only weakly coupled, and the syn- 
chrotron ages for radio-filled cavities are generally not 
good estimates of the dynamical ages (although in very 
young sources with ages much less than those in our sam- 
ple, the spectral ages are consistent with the dynamical 



spectral ages are roughly consistent with the buoyancy 
or sound-speed ages. 

4.6. Magnetic Field and Particle Content 

Because the magnetic field and particles are intimately 
connected through the observed synchrotron radiation, 
we discuss them together in the following sections. Ob- 
servational evidence for magnetic fields outside of the 
radio sources, but within the hot thermal gas in galaxy 
clusters, comes from Faraday rotation measurements of 



X-ray data 1 


e.g., Dreher et al. 1987 


Owen et al. 


1990 


Burns et al. 


1992 Ge & Owen 1993 


Taylor k \ 


J erley 



Radiative Efficiency and Content of Radio Sources 



9 



f993 


Ge & Owen|f994; 


Taylor et al.||1994 Clarke et al. 
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Taylor et al. 2002 


1. These studies show that the 



central cooling how regions have magnetic field strengths 
of 5 to 10 /xG, with a coher ence scale between 5 and 10 



for a given value of fc, the magnetic field required for 
pressure equilibrium, B p , can be determined from the 
X-ray and radio properties of a lobe. 
Considered as a function of B, the right hand side of 



kpc (jCarilli fc Taylor|[2l302j ) . Current theory posits that equation [H] has a minimum (for B = 3~ 2/7 B cq , equation 
AGN are responsible tor the magnetization of the ICM, F7I). so that there is a minimum Dressure consistent with 



AGN are responsible for the magnetization of the ICM, 
and they may also have an imp act on large-scale struc- 
ture and galaxy formation flKronberg et al. 112001 Carilli 
fc Taylor||2002| |Benson et~l.||2003| ). 

Because X-ray and radio data probe different proper- 
ties of the lobes, they can be used together to increase our 
understanding of the lobe contents. Equation [T] governs 
the total energy content of the lobes, comprised of mag- 
netic field and particle contributions. However, before 
the discovery of X-ray cavities, the total energy (the left- 
hand s ide of the equation ) could only be estimated indi- 
rectly (Owen et al. 2000). Cavities provide direct mea- 
surement of the total me chanical energy to within factor 
of a few (see Section 4.7 1. The other factors in equation 
[l] are constrained either by X-ray or radio data, with the 
exceptions of fc and B. By further assuming pressure 
balance and equipartition, we can constrain these two 
remaining parameters. 

Other assumptions, e.g., of source uniformity, are re- 
quired to fully determine the relationships between in- 
trinsic and observed properties of radio lobes. The ad- 
equacy of these assumptions is not evident in advance, 
making it is unclear how well the resulting relationships 
reflect reality. Observed properties can also be used in a 
variety of combinations to constrain intrinsic lobe prop- 
erties. In the following sections we consider a range of 
approaches to constraining lobe contents from their ob- 
served properties, generally relaxing assumptions as we 
proceed. 

4.6.1. Do All Lobes Have the Same Composition? 

Properties of the radio lobes, particularly th eir pres- 
sures, can be constrained by X-ray data (e.g., Blanton 



et al.[2001 1. X-ray cavities are assumed to be near pres 
sure equilibrium with their surroundings, since they are 
observed to be long lived and often surro unded by cool 
rims which d o not appear to be shocked (McNamara & 
Nulsen|[2007l ). The pressure in a radio lobe is the sum of 



the particle pressure, pp, and the magnetic pressure, Pb- 
In equilibrium, the total pressure, pp+pb, should match 
that of the surrounding thermal gas, p t h- Here we assume 
Pb = B 2 /8ir, which is strictly valid only for expansion 
perpendicular to the magnetic field. More generally, the 
effective magnetic pressure depends on the structure of 
the magnetic field in a lobe. For example, if the field 
is tangled and isotropic, the effective mag netic pressure 
would be pb — B 2 /24tt (e.g., Hughes|l991 1. Uncertainty 
in this relationship is a source of uncertainty that affects 
estimates of B p in particular. 

With pb = B 2 /8tt, the lobe pressure is given in terms 
of the particle and field energies by 

(1 + k)E e E B . . 

Pp + PB = ^v^ + W (5) 

in the notation of equations [l]- [3] Using equations [2] and 
[3] to express Eb and E e in terms of the magnetic field 
strength, this equation determines lobe pressure in terms 
of B. Thermal pressures fo r the gas can be co mputed 
from the X-ray observations ( Rafferty et al.|2006 1. Thus, 



[7 ) , so thac there is a minimum pressure consistent with 
observations. For any thermal pressure larger than the 
minimum, equation [5] has two solutions for B. For the 
smaller solution, the particle pressure dominates and, for 
the larger, the magnetic pressure dominates. Since the 
smaller solution typically gives unreasonably small mag- 
netic field strengths (the majority of them smaller than 
1 /^G), we chose the solution for which the magnetic 
pressure dominates. Very low magnetic field strength 
can also imply excessive inverse Compton X-ray emis- 
sion from a lobe, violating con straints obtained from the 
gamma-ray background (e.g., Stawarz et al. 2006 find 
that B is unlikely to be smaller than 10 /zG on average 
in systems such as these) . Because the magnetic pressure 
obtained in this way is close to the thermal pressure, B p 
is relatively insensitive to k. Values of B p in Table [3] were 
calculated by this method, using the gas pressure at the 
radius of the center of a lobe, and assuming k = and 
= 

The widel y used the minimum ener gy condition (first 
proposed by Burbidgc 1956) requires flPacholczykfl970} 

E B = \{l + k)E e , (6) 

where k is the ratio of the energy of heavy particles 
to that in electrons, E e . Thus it implies approximate 
equipartition between the energy of the magnetic field 
and that in relativistic electrons and heavy particles. Un- 
der the minimum energy condition, the magnetic field 
strength is 



B, 



eq 



B ct [{l + k)/<p] 



2/7 



where 



B r 



(R \2/7 2/7, \ r 2/7 T/ -2/7 

(6rr) ' Cl < 2 {a,v 1 ,v 2 )W a A V ' . 



(7) 



(8) 



This minimum-energy magnetic field strength estimate 
depends on the radio luminosity integrated between two 
fixed frequencies (10 MHz and 10000 MHz). Analyses 
using a fixe d low energy cutoff for the electrons ( Brunetti 



et al.|[l9~9~7[) or kn owledge of the acceleration mechanism 
( |Beck fc Krausc 2005) suggest that low energy particles 



may make a greater contribution to the pressure than we 
find. However, throughout this paper we use the classical 
equipartition magnetic field strengths. 

We now consider whether all lobes can have the same 
composition (i.e. the same value of k), assuming that the 
minimum energy (equipartition) condition applies. 

The equipartition magnetic field strengths given in Ta- 
ble[3land Table[4](-Boq, B cq 327) are computed from equa- 
tion^] above using the volumes from X-ray observations 
and 327 MHz radio maps respectively, assuming no heavy 
particles (fc = 0) and a volume filling factor, (</>), of unity. 
The assumption of fc = is the limiting case, where the 
energ y in the heavy particles is negligible (see De Young 
2006]) . The choice of = 1 is justified by X-ray obser- 



vations of cavities and shocks, since the cavities must be 
mostly empty of thermal g otherwise they would 

not be evident in X-ray images QMcNamara et al.|[2000] 



Fabian et al.|2000||Blanton et al.|2T)01| ). This lack of ther- 
mal gas and the requirement of pressure support implies 
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TABLE 4 

Magnetic field, synchrotron age, and particle content using the lobe sizes from 327 

MHz radio observations. 



Name 


B cci 327 

( M G) 


t eq 3Q - 
(10° yr) 


1 + fc oq327 b 


1 + ^'bouy327 b 


^eq,p327 


A133 


8.1 ± 0.2 


> 6.8 


800 ± 


300 


310 ± 


20 


55 


± 


5 


A262 


3.9 ± 0.2 


11.3 ± 0.1 


3100 ± 


1600 


3600 ± 


1600 


38 


± 


6 


MS 0735.6+7421 


3.9 ± 0.1 


10.8 ± 0.1 


2600 ± 


800 


320 ± 


90 


36 


± 


4 


Hydra A-outcr 


3.1± 0.4 


25 ± 6 


26000 


2200 


57 


± 


6 


M84 


24 ± 1 


< 0.32 


NS' 




NS C 




30 


± 


10 


M87 


16 ± 1 


< 0.57 


42 ± 


4 


40 ± 


10 


49 


± 


4 


Centaums 


10 ± 3 


< 1.1 


230 ± 


10 


110 ± 


20 


50 


± 


10 


A2052 


8.6 ± 0.3 


3.7 ± 0.2 


350 ± 


40 


320 ± 


20 


46 


± 


2 


MKW 3S 


6.7 ± 0.1 


> 8.8 


900 ± 


100 


180 ± 


30 


46 


± 


3 


A2199 


10.0 ±0. 3 


> 4.7 


440 ± 


40 


230 ± 


20 


58 


± 


2 


Cygnus A 


54 ± 4 


0.15 ± 0.01 


4.2 ± 


0.2 


NS C 




86 


± 


6 


A4059 


6.3 ± 0.2 


9.4 ± 0.2 


800 ± 


100 


600 ± 


100 


43 


± 


2 



a Assumes k — and <j) — 1. 

b Assumes 0—1. 

c No physical solution. 

a high filling factor. Additionally, the energy represented 
by the observed shocks requires that the radio lo bes dis- 
place most of the vol ume that they occupy (see McNa- 
mara & Nulsen 2007 1 . Values in Table [3] are computed 
using cavity volumes from X-ray observations. Those in 
Tabled] use lobe volumes from the 327 MHz radio obser- 
vations. 

In Figure [3] left, we compare the equipartition mag- 
netic field strength, B cq , to that required for pressure 
equilibrium, B p , for k — and (j> — 1. For all of the ob- 
jects, B cq cannot support the cavities against the ICM 
pressure, although Centaurus, Cygnus A, and the inner 
lobes of Hydra A are within a factor of two of equilib- 
rium. These systems are three o f the five that were clas- 
sified as radio-active in Section 14.31 For the other two 
(M84 and M87), the thermal pressure is too low to give 
a solution for B p , so they are not plotted. For these ob- 
jects the minimum lobe pressure is too large and they are 
likely in a driving stage where the lobes are still being 
inflated (see the discussion in the next section) . Whereas 
the magnetic field required for pressure equilibrium, B p , 
is insensitive to k as discussed above, the equipartition 
magnetic field does depend on k. For example, the esti- 
mates of B cq would increase by factors of ~ 4 and ~ 7, 
if 1 + k = 100 and 1 + k = 1000, respectively. Assum- 
ing 1 + k = 100 (Figure [3] right), many objects shift 
closer to equality. However, we note that for this value 
of k there are many systems absent from the figure, since 
the pressure equilibrium equation does not have a solu- 
tion. If 1 + k — 100 in the latter systems, they could not 
produce the observed synchrotron radiation unless their 
pressures exceeded the surrounding gas pressure. 

There is no single value of k for which the equipartition 
magnetic field strengths are equal to the field strengths 
required for pressure equilibrium across our entire sam- 
ple. Therefore, either the lobes do not have the same 
composition or equipartition does not apply. We investi- 
gate these possibilities further in the following sections. 

4.6.2. Varying Lobe Composition 

If the composition of the radio lobes, i.e., k, is vari- 
able, and the lobes are close to equipartition, then k 
can be determined by applying the equipartition and 
pressure equilibrium conditions together. By requiring 
the equipartition magnetic field (equation [7]) to be equal 
to the magnetic field required for pressure equilibrium 



(equation |5| , we can determine k as 



V 



6-KL lad ci 2 {a,v u v 2 ) V 13 



72rr 



-Pth 



7/4 



(9) 



In applying equation [3] which gives the ratio of electron 
energy to synchrotron luminosity, it is assumed that the 
radio observations measure the power do wn to the low 
frequency cutoff, which is not always true ( |Hughes|1991 1. 
For a > 1/2, the bulk of the electron energy resides in low 
energy electrons, so that, if the low frequency cutoff is 
unknown, the factor ci 2 (a, v 2 ) in equationplis poorly 
determined. This can contribute an order of magnitude 
or more of uncertainty to values of k given by equation 
[9j By contrast, under the equipartition and pressure 
equilibrium assumptions, the magnetic pressure pb = 
(9/13)pth and the total particle pressure p p = (4/13)p t h, 
so that they and the magnetic field are unaffected by this 
uncertainty. 

We have computed the heavy particle content, k cq , us- 
ing equation M (i.e., the heavy particle content of the 
lobes requiredto match their pressure to the surround- 
ing gas pressure under equipartition). Values of fc oq in 
Table [3] are computed using cavity volumes from X-ray 
observations. Values of fc eq 327 in Table H] use 327 MHz 
lobe volumes from radio observations. The filling factor 
4> = 1 in all cases. The gas pressure is measured at the 
radius of the center of a lobe. Under these assumptions, 
we find values of k from ~ 1 to 4000. 

From theoretical estimates k is predicted to range be- 
tween 1 and 2000, dependi ng on the mech anism that is 
generating the electrons ( Pacholczyk|1970 1. If these pre- 



dictions are correct, the high values of k eq determined 
here imply either that equipartition does not apply in 
the lobes of some of these radio sources or that there 
is additional pressure support (e.g., from hot thermal 
gas). The latter scenario needs to be investigated f urther 
with deep Chandra images (see e.g., |Gitti et al.j2007) ?). 
Alternatively, there may be substantial variations from 
source to source in the low energy cutoff of the electron 
energy distribution (i.e., in c\ 2 (a, v\, v 2 )). The value of 
k might also be affected by other differences in intrinsic 
radio source history, e.g., by synchrotron aging. 

The magnetic field strength corresponding to k oq , 
B cqp , is determined from the equipartition requirement 
Pb = (9/13)pth, with p t h measured at the radius of 
the center of the X-ray cavity (see Table [3l . Values of 
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Fig. 3. — Left: Equipartition magnetic field strength vs. the field strength required for pressure equilibrium (fc = 0, <f> ■ 
line denotes equality. Right: The same, but for 1 + k = 100. Symbols are the same as in Figure [2] 

■B eq ,p327 (Table [3} differ only because the gas pressure is 
different at the cen ter o f the typically larger 327 MHz 
radio lobes (Section |4.8| . In general, the magnetic field 
strengths determined from the radio lobe volumes are a 
little smaller than those determined for the X-ray cavi- 
ties, as the pressures for the former case are determined 
at slightly larger distances from the cluster centers. The 
equipartition magnetic field strengths are generally a lit- 
tle smaller than those required for pressure equilibrium 
(Bp) and both are larger than the clus ter-wide magnetic 
field estimates ( |Carilli fc Taylor||2002 1 . 

Our range of estimates tor the ratio of heavy-particle 
energy to electron energ y ~ 1 - thousands, is similar to 
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1). The dashed 
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that of Dunn & Fabian (2004) and iDunn et al. (20051 



However, for some systems our numbers differ by orders 
of magnitude. There are several reasons for the discrep- 
ancies. Because k is inversely proportional to the total 
radio luminosity, the differences can be due in part to 
differe nces in the total radio lu minosity. |Dunn fc Fabian| 



• * 


1 1 1 1 

f 




1 : 










<-# 










• 




















A 



























100 



1000 

v c (MHz) 



10000 



( 2004[ ) and |Dunn et al. ( |2005[ ) used the spectral indices 
from the literature and the 5 GHz radio flux in order 
to calculate the total radio luminosity. In contrast, we 
use measurements of the lobe fluxes at multiple frequen- 
cies to determine the spectrum and the total luminos- 
ity, which is generally more accurate and less sensitive 
to aging affects. Furthermore, we note that for some 
systems the c avity sizes and positions used by |Dunn fc 
Fabian| J2004| and |Dunn et "al~1 fl2005|) differ from ours 
qUirzan et aT]2 004 Ra flerty et al .|200b|). The differences 



in sizes and positions will further affec t the volume and 
buoyancy-age calculation (Section 4.4 1, both important 
in the derivation of k. 

The large range of k may be due to the range in in- 
trinsic radio properties such as the synchrotron age. To 
investigate this possibility further, in Figure H] we plot 
fc eq versus the break frequency. We notice a dichotomy 
between high and low break frequency systems, such that 
objects with high k have a low break frequency, as might 
be expected if the low-break- frequency objects lost their 
high-energy electrons, increasing the ratio between heavy 
(total) particle energy and electron energy over time. Be- 
cause no clear trend is present, aging is unlikely to fully 



Fig. 4. — The ratio of heavy-particle energy to electron energy 
for which the equipartition magnetic field strength equals the field 
required for pressure equilibrium versus the break frequency. Sym- 
bols are the same as in Figure [2] 



explain the range in k that we find. 

However, since the bulk of the energy resides in low en- 
ergy electrons (for a > 0.5), the total energy of the elec- 
trons can only change significantly when electrons near 
to the low energy cutoff begin to age significantly. Until 
this stage, k will be little affected by synchrotron ag- 
ing. When the lowest energy electrons have aged signifi- 
cantly, the break frequency should have evolved to below 
the low frequency cutoff. Since we measure a break fre- 
quency above the cutoff for the majority of our sources, 
either k should be largely independent of the age, or the 
break frequency does not directly reflect the age. Also 
we note that k is low in the radio-filled sources. If we 
assumed that k starts close to unity and evolves due to 
synchrotron losses alone, then 90% of the remaining elec- 
tron energy would have to be lost to boost k by an or- 
der of magnitude or more. That would require every 
source with large k to be older than the synchrotron loss 
time at the low energy cutoff. Nor can adiabatic expan- 
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sion explain the high values of k. If the heavy particles 
are relativistic, adiabatic expansion does not change k, 
and if they are non-relativistic, k is reduced by expan- 
sion, as adiabatic expansion favors relativistic particles 
(E p /E e cx V-V3). 

The trend in Figure [4] and the large values of k sup- 
port the idea of a h eavy jet (for a co mplementary view 
on this subject, see |De You ng 2006) or entrainment of 
heavy particles (see e.g., |Rossi et al. 2004| |Dunn et al.| 
2006a). Entrainment may play a role by redistributing 



jet energy from electrons to protons or suprathcrmal par- 
ticles, increasing k and decreasing radio luminosity as the 
electrons lose energy. It may explain how objects with 
the same break frequency can have different k values. 

Lastly, by equating the synchrotron age deriv ed fr om 
the lobe's radio spectrum (equation E] in Section |4~4t to 
the X -ray-derived age for the cavity (from Rafferty et al. 
2006 1 , we can place further limits on the magnetic field 



strength and lobe content. The resulting relation deter- 
mines the magnetic field strength in a lobe, -Bbuoy (Table 

t, required to make the synchtron age consistent with 
e age estimated from X-ray observations (the solution 
for Sbuoy with B > ^/^IB^fi is used). Note that, at least 
in some cases, the buoyancy age is likely to be an overes- 
timate of the ag e of a cavity ( |Wise et aL 2007 McNamara 
fc Nulsen|[2007| , causing Bbuoy to be underestimated. 

Also assuming pressure equilibrium, so that the total 
pressure of a lobe (equation pi) equals the surrounding 
gas pressure, pth, taking B = i?buoy then gives 



buoy 



r>2 

^buoy 
87T 



3VB, 



3/2 
buoy 



LrndCl2 



(10) 



This approach, which combines the synchrotron age with 
constraints from X-ray observations, determines k with- 
out the n eed to assume equipartition, cf. Sections |4.6.1| 
and 4.6.2 Values of Bbuoy (Table [3]) are used in equation 
10 to obtain fcbuoy for the X-ray cavities (Table k3| and 
^buoy327, using the 327 MHz lobe volumes (TableTO. 

We note that for some lobes and some magnetic field 
estimates in Tables [3] and [4] the magnetic pressure of the 
lobe exceeds the thermal pressure of the surrounding gas, 
allowing no physical solution for k. For -Bbuoy, this may 
be because t^ uoy overestimates the age of the cavity (due 
to projection or because the cavity is driven). In other 
cases, the magnetic pressures may simply be overesti- 
mated (such that pb > (9/13)p t /j) within the large un- 
certainties in the equipartition assumption (discussed at 
the beginning of this section). Given these qualifications, 
if all the assumptions that go into these calculations are 
correct, lobes for which the thermal pressure is too low 
to allow pressure equilibrium mus t be over-pressured and 



therefore in a driving stage (e.g., Owen et aL"||2000 1 



4.7. Scaling Relations between Synchrotron Power and 
Jet Power 

While deep X-ray observations are critical for measur- 
ing the AGN output (cavity and shock energy), these 
observations are difficult to obtain in many clusters, par- 
ticularly at higher redshift. However, it is important 
to quantify the history of AGN output in clusters for a 
better understanding of the feedback process in galaxy 
formation and cluster preheating (e.g., Magliocchetti & 



Briiggcn 2007; Bes t et aL||2007 l. This goal requires sen- 
sitive measurements of both the current rate of energy 
injection and the history of injection with cosmic time. 
Radio observations, properly calibrated to the total en- 
ergy, may be a useful proxy for X-ray observations as 
a tracer of the energy injection. Our goal in this sec- 
tion is to use the observed cavity and radio properties to 

calibrate the radio-to-total jet power. 

The first attempt at this was presented in |Brrzan et al.| 

l /2 

(2004), where we found that B r adio ~ -Pcav , but with 
large scatter. This scatter and the unknown temporal 
behavior of the radio emission limits the usefulness of 
this relation to ac curat ely probe jet power. Nevertheless , 
Best et al.| (|2007[) and |Magliocchetti fc Briiggen| ([2007 1 



adopted the Birzan relation to study feedback in elliptical 
galax ies using 1.4 GHz ra dio power alone. Like the stud- 
ies o fjBirzan et al.| ( |2004| and |Rafferty et al.| ( [2006| , they 
found that AGJN were able to quench cooling in many sys- 
tems. However, an additional source of heat is implied 
in massive systems. This conclusion is uncertain because 
it takes no account of the large scatter and its origin in 
iradio _ -Pcav relation. For example, one powerful outburst 
in any system can dominate the AGN heat input over its 
entire lifetime. In this case, the variance in the relation 
would be as important as the mean. Furthermore, the 
slope and scatter depend on frequency. Frequencies be- 
low 1 GHz are most sensitive to feedback power. A more 
accurate relation that can be applied to data from all-sky 
surveys would p rovide the basis for understanding radio- 
mode feedback ( |Croton et al. 2006| |Sijacki fc Springe!] 



2006 1 and its impact on galaxy formation and evolution. 
Additionally, accurate estimates of the radiative efficien- 
cies are important for jet and accretion models. With 
these goals in mind, in this section we use our sample of 
cavity systems to examine the relationship between the 
radio power and the cavity power. 

In Figure [5] we plot the jet power versus the bolo- 
mctric radio luminosity for the total source (left) and 
the lobes (right) of the radio sources. This figure shows 
the same trend between bolometric radio l uminosity and 
cavity power as seen in jBirzan et al.| ( |2004[ ) . The most 
radio-luminous objects generally have the largest cavity 
power. This trend is shared by both radio-filled cavities 
and radio-ghost cavities, but the ghost cavities tend to 
have a higher ratio of jet power to radiative power. 

To quantify the trend, we perform an ordinary least 
squares regression to fit a linear function to the loga- 
rithms of the data. The best-fit line for the total source 
luminosity is given by 

logP cav = (0.48 ± 0.07) logi radio + (2.32 ±0.09), (11) 

with a scatter (standard deviation) of ss 0.83 (see Figure 
[5j left). The best-fit line for the lobes only (Figure [5j 
right) is 

logP cav = (0.59 ±0.08) log L radio + (2.11 ±0.11), (12) 

with a scatter of f=a 0.65 dex. In these relations, both the 
cavity power and the bolometric radio luminosity are in 
units of 10 42 erg s _1 . The relation between the jet power 
an d total radio power is similar to our previous finding 
in Birzan et al. ( |2004[ ), where we used inhomogeneous 
literature data in order to quantify the radio power. 

When we use only the five radio- filled sources, those 
with break frequencies above 1400 MHz, the best fit line 
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Fig. 5. — Cavity (jet) power versus bolometric radio power for the total source (left) and lobes only (right). The symbols and wide error 
bars denote the values of the cavity power calculated using the buoyancy timescale. The short and medium-width error bars denote upper 
and lower limits of the cavity power calculated using the sound speed and refill timescales, respectively. Each point represents the sum of 
cavity powers of each cavity type. Filled symbols denote radio-filled cavities, and open symbols denote radio-ghost cavities. Dotted lines 
denote diff eren t r atios of cavity power to bolometric radio luminosity. The dashed line show s th e be st-f it power law for the entire sample 
(equations |11| and |12| ; the solid line shows the fit for the radio- filled systems only (equations |13| and |14) . 

in units of 10 erg s _1 is given by 



log P cav = (0.70 ±0.12) log L radio + (1.40 ±0.17) (13) 

for the total source (the solid line in Figure [5j left) and 
for the lobes only (the solid line in Figure [5] right) by 

logP cav = (0.69 ±0.12) log L radio + (1.46 ±0.17). (14) 

For the radio-filled sources, the jet power scales tightly 
with both the total and lobe radio power, with a scatter 
in both cases of only 0.31 dex. 

For comparison with monochromatic radio data, we 
also plot the cavity power against the 1400 MHz and 
327 MHz powers of the total source in Figure [6j The 
monochromatic radio power was calculated as P„ = 
AttD\S„(1 + z) a ~ 1 , where we have assumed a power-law 
spectrum (S u oc v~ a ). The best-fit relations are (with 
P cav in units of 10 42 erg s _1 and P327 and P1400 in units 
of 10 24 WHz- 1 ): 

log Pcav - (0.62 ± 0.08) logP 327 + (1.11 ±0.17) (15) 

log P cav = (0.35 ± 0.07) logP 140 o + (1.85 ± 0.10). (16) 

The scatter about the 1400 MHz relation for the whole 
sample is 0.85 dex, compared to 0.69 dex for the 327 
MHz relation (an improvement of « 20%). In terms of 
the scatter, the monochromatic 327 MHz power is al- 
most as accurate for inferring the c avit y power as the 
bolometric lobe luminosity (equation 12 1. We note that 



the ratio of cavity power to radio power ranges between 
a few and a few thousand. About half of our systems 
have ratios in the range 10-100, with the rest being on 
average vastly larger (radiative efficiencies of < 0.001). 
The median value for the ratio of cavity power (ApV/t) 
to radio power is 420 when we consider the total radio 
luminosity, and 141 when we consider only the lobe ra- 
dio luminosity. However, the mean values for the ratio 
of cavity power to radio power are 4200 (total radio lu- 
minosity) and 595 (lobe radio luminosity). The scatter 
that we see in Figure [5] is much larger than the range in 
theoreti cal estimates, which generally pred ict a ratio of 



10-100 fDe Young|1993||Bicknell et al.|1997[ ). Asaresult, 
the total or lobe radio luminosity is not a good predictor 
of the total jet power. It is important to note that we 
did not include the contribution of shocks in the cavity 
power calculation for any of the objects. Deep Chandra 
X-ray images show that the lobes create mild shocks dur- 
ing their exp ansion, with a Mach number lying between 
1.2 and 1.7 ([Fabian et al.||2006| IMcNamara et al.||2005 



.\\^KjKiuaiu\ |Forman et al.||2005| | Wise et a 



Nuls en et al.||2005a|b[ . 

20071 IForman et al.||2007| IMcNamara & Nulsen||2007 1 



the slopes of the best-fit relations differby almost a fac- 
tor of 2, whereas one might expect them to have similar 
slopes (if P327 scales roughly with Pi 400)- However, we 
are not attempting to uncover the true, underlying rela- 
tion between cavity power and radio luminosity. Rather, 
we wish to find the best relation to use for predictive pur- 
poses. Therefore, the relation that results from a simple 
ordinary least square s regression is the most appropriate 
(for a discussion, see Isobe et al.||1990[ ). 

The lar ge scatter about the se trends is similar to that 
found in Birzan et al. (2004). In Figure [5] we see that 



In cases where shocks are present, the cavity powers are 
lower limits to the jet powers, and as a consequence the 
radiative efficiencies are overestimated. 

Radio-source aging may be responsible for some of the 
scatter in this plot: in general, the most radiatively ef- 
ficient systems (e.g., Cygnus A and M87) are also the 
youngest ones and may still be driven by the radio source. 
Conversely, the most radiatively inefficient systems are 
generally the older ones. In this scenario, the sources 
start with high radiative efficiencies when they are young 
and, as they age, their efficiencies decrease. However, 
since the radiative efficiencies depend mainly on the lower 
energy electrons which have a very long radiative life- 
time, the radiative efficiencies of the sources in our sam- 
ple should not be strongly affected by aging. 
To further investigate whether the scatter is due to ra- 
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dio aging, in the left panel of Figure [7] we plot the ratio 
of cavity power to radiative power (the inverse of radia- 
tive efficiency) versus the break frequency. Only a few 
objects are radiatively efficient (e.g, M87, Cygnus A, and 
Hydra A - inner). If the scatter in Figure pj is a radio- 
aging effect, one would expect to see a correlation be- 
tween radiative efficiency and the break frequency, such 
that those objects with high efficiencies have high break 
frequencies (and are therefore younger). Figure m left, 
shows a segregation between young (radio-filled) objects, 
which may be still in a driving stage (e.g., M87, Cygnus 
A, etc.), and the older (ghost) objects (e.g., Hydra A 
- outer, A1835, etc.) and a tendency for objects with 
higher break frequencies to be radiatively efficient. This 
segregation is similar to the one in Figure [4] This result 
is not surprising since k is approximately proportional to 
the ratio between cavity power squared and the bolomet- 



ric radio luminosity (see equation [9]) . As a result, the k 
values are indicators of the radiative efficiencies for the 
radio sources. As we discussed in section |4.6| the scatter 
in the radiative efficiencies may be due to a combina- 
tion of aging and entrainment, which will increase k and 
reduce the bolometric radio luminosity. 

A similar segregation is seen in the right panel of Figure 
[7J where we plot the ratio of cavity energy to radiative 
power versus the break frequency. The ratio of cavity en- 
ergy to radiative power represents the timescale for the 
AGN to radiate away its cavity energy via synchrotron 
radiation. Since many of the objects have a very long 
timescale, longer than the Hubble time, the synchrotron 
radiation is generally a negligible fraction of the total 
energy budget currently. Only for a few objects is radi- 
ation important (e.g, M87, Cygnus A, and M84). As a 
result, aging does not appear to be the only factor that 
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contributes to the large scatter in the radiative-efficiency 
plot. 

However, the scatter may be reduced if the dependence 
of the radiative efficiencies on break frequency is taken 
into account. To this end, we fit the equation logP cav = 
A log L ra dio + B log vq + C to all objects with estimates of 
the lobe break frequency. For the purposes of this fit we 
treat the upper and lower limits on the break frequency 
as detections. The resulting best-fit relation is: 

log P cav = (0.53 ± 0.09) log L radio - 

(0.74 ± 0.26) log v c + (2.12 ± 0.19), (17) 

where the cavity power and lobe radio luminosity are in 
units of 10 42 erg s" 1 , and vq is in units of GHz. 



In 



versus 



Figure [8] we plot the cavity power 

(^radio) ' 53 (^c) -0 ' 74 - Figure [ijshows that accounting for 
the synchrotron aging greatly reduces the scatter about 
the best-fit line, from 0.64 dex (equation 12) to 0.33 dex. 
Therefore, knowledge of the lobe break frequency im- 
proves the scatter by ~ 50%, significantly increasing the 
accuracy by which one can estimate the jet power of the 
AGN when cavity data are unavailable. A relation with 



similar scatter was identified by |Merloni fc Heinz ( 2007 1 
using the 5 GHz core radio luminosity only, after correc- 
tion for relativistic beaming. 

Another factor that may contribute to the scatter in 
the radiative efficiency plot is the interplay between the 
magnetic field and the particle content. From Figure 
[31 we concluded that for k = 0, only in a few sources 
(the youngest ones) can the equipartition fields achieve 
pressure equilibrium. For the majority of the objects, 
in order for the equipartition fields to achieve pressure 
equili brium, a much larger v alue of k is required (Section 



4.6.11. However, in Section 4.6 we concluded that the 



large range in k is not due only to variations in magnetic 
field strengths, but may be due to differences in intrinsic 
radio properties. Because of the interrelations between 
magnetic field, particle content, and intrinsic radio prop- 
erties such as radio luminosity, the scatter in Figure [5] is 
not due exclusively to one of these factors, but may be a 



4.8. Radio Observations As a Tracer of Cavity Size 

Low-frequency radio observations appear t o be a good 
trace r of cavity activity over 10 s yr timescale s (Lane et al 
" Fabian et al.|| 2001[ ) 



2004 



Wise et al. 



2007 



For exam- 
ple, 327 MHz VLA radio maps have proved to be crucial 
for detect ing the outer, fainter cavities in sy stems like 



Hydra A (Wise et al. 
combination of a brig 



2007| |Lane et al Jj 2004), due to a 
iter source at 327 MHz (particu- 
larly important for steep spectrum sources, such as MS 
0735.6+7421) and the generally higher sensitivity of the 
VLA at low frequencies to larger-scale structure (though 
depending on the array this may not always be the case) . 
We note, however, that the system temperature of the 
VLA is higher at lower frequencies, so there is some trade 
off between the above effects and VLA's sensitivity to low 
surface brightness emission. 

In order to evaluate whether 327 MHz emission is gen- 
erally a better predictor of cavity size than 1400 MHz 
emission, we compare pV cav , where V cav is the volume of 
the cavities from the X-ray maps against pT^i4oo, where 
V1400 is the volume of the lobes from 1400 MHz radio 
maps, in Figure [9j left. We also plot pV cav against pV 32 7, 
where V327 is the volume of the lobes from 327 MHz ra- 
dio maps (Figure [9j right). The pressure on each axis 
is measured at the position of the radio lobe's center or 
cavity's center, and can therefore differ from one axis to 
another. In Table [5] we list the the lobe's size and the 
pressure at the lobe's center for both 327 MHz and 1400 
MHz radio lob e s. Th e X-ray pressures are taken from 
Rafferty et al. (20061. We include only those systems 



with resolved emission in this analysis. 

Figure [9] suggests that 1400 MHz emission in our im- 
ages is a good tracer of cavity size in most cases, but there 
is a significant scatter, with many points falling below the 
equality line. Assuming that P C av traces the minimum 
jet power, we note that the power we estimate from the 
1400 MHz emission underestimates the cavity power in 
some systems. On the other hand, at 327 MHz, all of the 
points are above the equality line, despite the fact that 
the radio data have a variety of noise levels and differing 
uv coverage and antenna spacings (array configurations), 
all of which affect the sensitivity of the resulting images 
to the lobe structure. This effect is likely due to steep- 
ening of the spectrum towards the edges of the lobes, a 
common property of many low-power radio sources (e.g., 
Klein et al.||1995| ). It therefore appears that, in general, 
cavity energies inferred from lower-frequency radio maps 
are less likely to underestimate the X-ray-inferred ener- 
gies. 

A number of objects lie above the equality line at both 
frequencies (e.g., Abell 4059, Perseus, Abell 2052, RBS 
797, and 2A 0335+096). In these cases, the radio lobes 
are larger than the cavities measured in X-rays. It is pos- 
sible that X-ray measurements have underestimated the 
sizes of these cavities, or that the radio plasma has dif- 
fused beyond the cavities. The latter explanation seems 
to be favoured for Abell 4059 if we think of it as an inter- 
mediate case in wh ich the X-ray cavities are fading away 
(iHeinz et al.l 



In Hydra A, the inner cavities that we measured from 
X-ray images are filled wit h 4500 MHz radio emission. 
However, from Wise et al. (20071 we know that Hydra 
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TABLE 5 
Radio lobe sizes and pressures. 



Name 


Component 


0,327 X b327 a 

(arcscc X arcscc) 


P327 b 
(10 -11 erg cm~ 3 ) 


ai.4 X 
(arcsec X 


fcl.4 C 

arcscc) 


Pl.4 d 

(10 11 erg cm 


A133 


total 


24 X 23 


17 ± 3 


31 x 


10 


10.7 ± 0.4 


A262 


cast lobe 


19 x 11 


9 ± 2 


16 x 


9 


10 ± 2 




west lobe 


26 x 11 


8 ± 2 


12 x 


8 


12 ± 2 


Perseus 


north lobe 


30 x 23 


100 ± 20 


20 x 


17 


120 ± 20 




south lobe 


35 x 33 


70 ± 20 


15 x 


14 


160 ± 20 


2A 0335+096 


total 


53 x 16 


24.7 ± 0.4 


59 x 


35 


23.5 ± 0.4 


A478 


total 


< 6.1 x 4.6 


130 ± 20 










north lobe 






1.6 x 


1.3 


170 ± 20 




south lobe 






1.9 x 


1.6 


170 ± 20 


MS 0735.6+7421 


north lobe 


46 x 30 


7.0 ± 0.8 


35 x 


23 


9.3 ± 1.1 




south lobe 


41 x 30 


8.2 ± 1.2 


38 x 


27 


8.7 ± 1.1 


Hydra A-outer c 


north lobe 


120 x 90 


2 f 










south lobe 


120 x 45 


5 f 








Hydra A-inncr g 


north lobe 


45 x 30 


23.0 ± 0.4 


34 x 


21 


23.0 ± 0.4 


south lobe 


45 x 30 


23.0 ± 0.4 


35 x 


20 


22.0 ± 0.4 


RBS 797 


total 


25 x 17 


63 ± 7 










cast lobe 






4.2 x 


3.7 


92 ± 7 




west lobe 






3.9 x 


3.2 


92 ± 7 


Zw 2701 


total 












Zw 3146 


total 












M84 


north lobe 


40 x 31 


2.7 ± 1.5 


39 x 


33 


2.4 ± 0.2 




south lobe 


42 x 26 


2.6 ± 0.2 


46 x 


26 


1.6 ± 0.2 


M87 


cast lobe 


23 x 20 


38 ± 4 


25 X 


21 


38 ± 4 




west lobe 


24 x 19 


45 ± 3 


22 X 


17 


45 ± 3 


Ccntaurus 


total 


72 x 43 


13.3 ± 0.3 










cast lobe 






10 x 


10 


27 ± 3 




west lobe 






13 x 


13 


24 ± 3 


HCG 62 


south lobe 






< 5.6 


< 2.8 


9 ± 1 


A1795 


north lobe 






3.4 x 


6.2 


60 ± 10 




south lobe 






4.7 x 


2.6 


60 ± 10 


A1835 


total 












MACS J1423. 8+2404 


total 












A2052 


total 


37 x 21 


12.2 ± 0.8 


32 x 


22 


11.1 ± 0.7 


MKW 3S 


total 


59 x 54 


12.2 ± 0.9 


30 x 


21 


11.1 ± 1.1 


A2199 


cast lobe 


31 x 17 


19.5 ± 1.1 


13 x 


4 


27 ± 7 




west lobe 


30 x 20 


19.7 ± 1.1 


9.0 x 


7.3 


27 ± 7 


Cygnus A 


cast lobe 


40 x 17 


36.9 ± 0.9 


32 x 


15 


33 ± 2 


west lobe 


39 x 17 


38.8 ± 0.9 


34 x 


17 


31 ± 2 


Sersic 159/03 


total 


52 x 14 


18.6 ± 1.3 








A2597 


cast lobe 


19 x 19 


36 ± 5 










west lobe 


25 x 17 


35 ± 6 










total 






10 x 


8 


40 ± 4 


A4059 


north lobe 


27 x 24 


11 ± 1 


23 X 


18 


9 ± 1 




south lobe 


35 x 28 


11 ± 1 


25 X 


20 


8 ± 1 



Projected semimajor and semiminor axes of th e lobes from the 327 MHz radio map. 
The X-ray pressure (from 



Rafferty et 



al.|2006| l at the center of the lobes measured from 327 MHz radio maps. 
ijees of th e lobes from the 1400 MHz radio map. 



radi< 



c Projected semimajor and semiminor axes i. . .... . ,.. 

The X-ray pressure (from [Rafferty et al. |2006[ l at the center of the lobes measured from 1400 MHz : 
c The outer cavities; the cavities called E and F in |Wise et al.| | |2007} . 
The pressure is from |Wise et al.N2007l. 

■ cavities measured by [BTrzan~et al.| < |2"004[ l and |Rafferty et al.| < f2~006| l ; called A and B in |Wise et al~| ( |2007t . 
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A has outer caviti es that are filled with 327 MHz ra- 
dio emission (from Lane et al. 2004 1. For Hydra A, we 



calculated two sets o f data: one for the i nner cavities 
using the pV cav from |Rafferty et al.| p006| and o ne for 



the outer cav ities using the X-ray information from | Wise] 
et al. ( |2007[ ). There are other objects similar to Hydra 
A where the energy inferred from the 327 MHz radio 
emission is well above the energy we measure from X-ray 
images (e.g., RBS 797). In these objects we expect to 
find outer cavities in deeper Chandra images which will 
increase their total pV cav estimates. 

5. CONCLUSIONS 

In summary, using an X-ray and radio study of 24 sys- 
tems with X-ray cavities, we found that these systems 
are in general poor radiators with a ratio of cavity (jet) 
power to radio power ranging from a few to a few thou- 
sands. We conclude that this scatter is probably due 
to several effects including aging and entrainment plus 
a combination of intrinsic differences in magnetic field 
strength and particle content. Accounting for the de- 
pendence of the radiative efficiencies on aging reduces 
the scatter (standard deviation) of the relationship be- 
tween radio luminosity and average jet power by ss 50%, 
improving the accuracy with which one can determine 
the total jet power of the AGN using radio observations 
alone. However, age estimates are costly to obtain, as 
an accurate estimate of the break frequency requires ob- 
servations at several different frequencies. When data 
at only a single frequency are available and in the ab- 
sence of an estimate of the break frequency, we found that 
low-frequency (327 MHz) monochromatic radio power al- 
lows a more accurate estimate of the jet power than the 
higher-frequency 1400 MHz powers (the scatter at 327 
MHz is w 20% lower than at 1400 MHz). Therefore, 
low-frequency radio observations may be the cheapest 
strategy for studying feedback in radio galaxies. 

Our data allow us to put limits on the particle con- 
tent and, through comparisons of the buoyancy and syn- 
chrotron ages, to test the assumptions that go into their 
determination. We found that the two age estimates are 
weakly coupled. We argue that the synchrotron ages un- 
derestimate the true ages for the younger systems (with 
high break frequencies), which likely are actively being 
inflated by the radio source and which generally have a 
lower value of k. On the other hand, for the older sys- 
tems (with low break frequencies), the cavity ages and 
synchrotron ages are closer to being in agreement, and k 
is generally much higher. 

Using both radio and X-ray data we placed limits on 
the magnetic field strengths and particle content in the 
lobes (k, defined as the ratio between the heavy-particle 
energy and the electron energy) . Assuming a lepton jet 
(k = 0), we found that the magnetic field strengths re- 
quired for pressure equilibrium arc in general much larger 



than both the equipartition magnetic field strengths and 
those derived by requiring the synchrotron age to be 
equal to the buoyancy age. Therefore, either equipar- 
tition does not apply in the lobes of the radio sources or 
there is additional pressure support from other sources. 
If we allow for heavy jets {k > 0), equating the equipar- 
tition magnetic field strength to that required for pres- 
sure equilibrium gives values of k between approximately 
unity and 4000 a ra nge similar to that found by |Dunn| 



fc Fabian (2004) and|Dunn etaL (20051. The range in k 



may be due to the effects of aging and entrainment 

Our multiwavelength radio data allow us to better clas- 
sify the cavities based on their radio emission. Motivated 
by the traditional classification system based on the pres- 
ence or absence of 1400 MHz emission in the cavity, we 
used the break frequency to classify the cavities as ghost 
{v c < 1400 MHz) and active (v c > 1400 MHz). Based 
on this classification, only M84, M87, Hydra A, Cygnus 
A, and Centaurus have radio-filled cavities. However, 
we note that since our systems have a continuous range 
in break frequencies, any separation of the cavities into 
discrete categories is somewhat arbitrary. 

Lastly, by comparing the cavity sizes inferred from 
the X-ray maps with those inferred from 1400 MHz and 
327 MHz radio maps, we found that the sizes from our 
1400 MHz images often under-predict that cavity energy, 
whereas the lobes sizes measured from our 327 MHz im- 
ages generally equal or exceed the X-ray cavity sizes. 
Where the low-frequency lobe emission appears to ex- 
tend far beyond the X-ray cavities, there may exist faint 
outer cavities which would only be visible in very deep 
X-ray images (such as the outermost cavities in Hydra 
A). 

The discovery of X-ray cavities in the hot atmosphere 
of clusters has made it possible to estimate the radiative 
efficiencies of the radio sources, which can be much lower 
than the theoretical estimates (as low as 0.0001). There- 
fore, most of the AGN's energy in such radio galaxies is 
dumped directly into the ICM through the X-ray bubbles 
and shocks, likely preventing the overcooling of gas below 
~ 2 keV which would lead to an over-abundance of ex- 



tremely bright galaxies at the center s of clusters (Sijacki 
fc Springel|[2006| |Croton et al.||2006[ ) 
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APPENDIX 
RADIO IMAGES 



Figure 10 presents the radio contours for each image listed in Table [7J Figure [TT] shows the resolved radio emission 



superimposed on the smoothed Chandra X-ray images. In this figure, systems are ordered by decreasing P C av/^radii 
(i.e., increasing radiative efficiency), and images for each system are ordered by increasing frequency. 

DETAILS OF RADIO OBSERVATIONS AND IMAGE PROPERTIES 

Table [6] lists the details of our observations for each object in the sample. Table [jj lists the basic radio image 
properties such as resolution, noise, position angle for the beam at each frequency, anathe arrays from which data 
were used to make the image. Also listed are the references to published images. 



Fig. 10. — Examples of the radio contour plots. See http://www.phy.ohiou.edu/~birzan for a complete version. 
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TABLE 6 
Summary of the radio observations. 





Frcqucncy a 


Bw 






Duration 






Name 


(MHz) 


(MHz) 


Array 


Date 


(hh:mm:ss) 


Flux calibrator 


Phase calibrator 


A133 


327.5/333 


3.125 


A 


30 November 1992 


00:68:00 


3C48 


0023-263 




328.5/321.56 


6.25 


AB 


9 June 2002 


02:19:53 


3C48 


0116-208 




1464.9/1385.1 


50 


A 


25 June 1998 


02:31:10 


3C286/3C48 


036-216 




1364.9/1435.1 


25 


C 


9 December 2002 


03:08:10 


3C48 


0116-208 




1364.9/1435.1 


25 


BC 


29 September 2002 


01:07:40 


3C48 


0116-208 




1364.9/1435.1 


25 


BC 


6/7 October 2001 


00:57:40 


3C48 


0116-208 




8085/8335 


50 


D 


3 August 1988 


00:38:20 


3C48 


2234+282 




8535/8785 


50 


D 


3 August 1988 


00:37:00 


3C48 


2234+282 


A262 


327.5/321.56* 


3.125 


A 


13 March 2006 


04:23:40 


3C48 


3C48 




1515.9/1365.1* 


25 


B 


14 May 2005 


02:13:10 


3C48 


0119+321 




8535.1/8485.1* 


50 


B 


14 May 2005 


01:38:50 


3C48 


0205+322 


Note. 


— For a complete 


version, 


sec http : //www . phy . ohiou . eduA 


^birzan. 







a Our observations arc marked with *. 



The integration time on source. 



TABLE 7 
Properties of the radio images. 





Frequency 




Resolution 


PA 


Rms noise 




Name 


(MHz) 


Array 


(arcscc X arcsec) 


(degree) 


(mjy/bcam) 


References 


A133 


330.25 


A 


10.3 X 6.6 


12.6 


1.72 






321.56 


AB 


34.4 x 18.5 


-23.1 


2.03 






1425 


A 


3.5 x 2.3 


65.1 


0.046 


21 




1400 


A, BC, C 


20.3 x 12.4 


-0.7 


0.060 


14 




8660 


D 


18.8 x 6.9 


-27.8 


0.033 




A262 


324.531 


A 


6.4 x 5.5 


-29.8 


0.55 






1365.1 


B 


4.3 x 3.5 


-48.6 


0.045 


17 




8460.1 


B 


0.7 x 0.6 


49.2 


0.021 




Perseus 


327.5 


A. B. C 


8.6 x 7.6 


-83.8 


7.6 


6 




1575 


A, B, C 


8.1 x 4.4 


86.5 


3.0 


6, 18 




8414.9 


BC 


2.0 x 1.2 


82.6 


15.0 




2A 0335+096 


324.531 


B 


21.7 x 16.4 


-37.4 


3.7 






1489.9 


C 


14.3 x 13.8 


-47.7 


0.050 


20 




4860.1 


D 


15.7 x 13.9 


6.7 


0.008 


20 




8439.9 


C 


4.7 x 4.7 


-45.0 


0.025 


20 



Note. — For a complete version, see http://www.phy.ohiou.edu/~birzan. 



